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I. SAMPLE FUNCTION AND DIAMETER MEASUREMENT 


1 . Sample Function Calculation 

The sample function of the PUDVM can be measured experimentally 
as a step velocity response. As the sample region scans a zero 
velocity/high velocity interface the recorded RMS power from the 
audio signal can be used to calculate the effective length of the 

sample function. For ^transducer ^^^vessel’ niethod provides 
an accurate measure of the sample function length since the experi- 
mental constraint Is that the step be abrupt compared to the resolution 
of the ultrasound system. For these determinations flow in the 
vessel is maintained at RE=4000 to insure a blunt profile and there- 
fore steep velocity gradient at the wall. The abruptness of the 
experimental test function can be estimated from the minimal velocity 
to which the instrument will respond and the velocity gradient dv/dr 
at the wall. Assuming a turbulent velocity profile: 

V = 
and 


dr % 


where is the peak velocity and r/R is the nondimensional radius, 
The velocity gradient at r/R equal to 1 is: 


dv 

d(ft 


r/R=l 


-4V 


The minimal sensible velocity of the current PUDVM design (McLeod 1974) 
is on the order of Icm/sec. A typical = 50 cm/sec. The abruptness 
of the step A(-^) = The 1/200 fractional radius step abruptness 
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is far steeper than the instrument resolution and can therefore be 
considered ideal. 

Figures 1 and 2 exhibit RMS voltage curves for a transducer 
diameter equal to the sample function is scanned 

across the vessel. Each successive curve is obtained by increasing 
the gate setting from 1 usec to 8 psec in 1 psec increments. The 
sample function is obtained as indicated by calculating the delay 
times of the rise of the RMS voltage from zero to its peak. To be 
most precise the voltage curves should be converted to power curves. 
For a gate setting of 1 psec the sample function is 1.8 psec implying 
a length Lg of 

Lq = -tq • 1500m/sec . s1n(60") , 

2 

Therefore gate settings can be easily found to provide sample function 
lengths of R/2, R, 3R/2, 2R for the experiments to be described. 

2. Measurement of Vessel Diameter 

Vessel diameter can be measured accurately by recording audio 
power during a sample function scan of the vessel cross-section. 

The half power points at the near and far slopes correspond to the 
locations of the walls of the vessel, since half power assumes the 
centroid of the sample function is at the wall -fluid interface. 

An example of a half for transducer R/2 for the four 

different gates is shown in Figure 3. Note that the measured diameter 
is 7.3 mm which is within S% of the actual vessel diameter. Note also 
the poor wall discrimination for a gate greater than R/2. 



4 ^ 15*0 



FIGURE 1 





RMS VOLTAGl’ FOR GATES 
• R/2, R, 3R/2, 2R 
i TRANSDUCER DIAMETER • R/2 
t RE ■ 4000 


FIGURE 3 
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FIGURE 7 
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3. Velocity Profiles 

Velocity profiles are obtained by manually scanning the sample 
function across the lumen and recording the output of the zero crosser 
(Histand, 1973). Typical velocity profiles for transducer R/2 for 
the 4 different gates are shown in Figure 4. 

4. Results of Power Scans 

Power scans in 1 psec increments across the vessel diameter 
were measured for the three transducer diameters; R/2, R, 2R, Data 
for transducer R is shown in Figure 5 and for transducer 2R in Figure 
6. The power curves for 2R were obtained with a flattened tube to 
provide an interface across the full transducer face. Figures 7 
and 8 show power scans for transducer R and 2R respectively for gates 
of R/2, R, 3R/2, 2R. 

5. Discussion 

For a narrow gate (<R) small transducer {<R) one can accurately 
^tiasure the vessel diameter. But as the gate widens from these values 
the power curves no longer plateau and it is impossible to estimate 
the half power points and thus the diameter. 
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II. VELOCITY PROFILES AND FLOW CALCULATION 


The purpose of these experiments was to assess the capability 
of the PUDVM as a velocity and flow measuring instrument in a simulated 
transmural application. The nondimensional parameters: transducer 
diameter, gate length, flow rate (Reynold's number), and vessel 
diameter were designated so that an investigator could judge the 
resolution of his measurements in an actual transcutaneous blood flow 
application. This experimenta’ data is to be compared with the 
theoretical data described on pages 6-46 (Daigle, 1974). The following 
parameter and protocol was carried out with the end goal to determine 
velocity and flow measurement accuracy by various data processing 
schemes. 

1 . Convolution/Bounda r y/Truncation Dialysis Tube - Transducer Experiment 

Parameters Controlled. 

1. 7.7 mm dialysis tube 200 DIA entrance length; normalize 
radius = 1 , for comparison of any vessel diameter, transducer diameter, 
PUDVM gate ratios. 

2. Steady Flow; Re=1000, Re=1500 (laminar), Re=3000 (turbulent). 

3. Transducer diameter (aluminum epoxy backing). 


Dia 

Normalized Diameter 

7.7 mm 

2R 

:3.85mm 

R 

1.98mm 

R/2 

Pulse length = 

8 cycles " x = 1 ysec 


5. Gate Settings (determined by power scans over flow/.no flow 


interfaces) . 
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length 

R/2 

R 

3R/2 

2R 

Data. 

1. RMS voltage curves (-power curves) for 1/2 pov/er point 
determination. 

2. Integrate profiles (standard computer program; diameter 
from max/mi n slope. 

3. Integrate truncated profiles to known and half power diameters. 

4. Small crystal wide gate calculation. 

5. Compare velocity values, particular attention to full 
illumination wide gate. 

Summary. 

(3 transducers )x (4 gates)x(3 Re values) = 36 runs 
Tabulated and plotted output. 
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2 . Di ameter Profile Integration Method 

A technique for measurement of blood flow involves determination 
of the blood velocity profile across a diameter of the vessel under 
study. Since this technique requires measurements of local blood 
velocity at discrete points across the vessel, the PUDVM is employed. 
The procedure usually followed with single gate instruments is as 
follows; The transducer is positioned at a known angle to the vessel 
and the sample volume located outside the near wall of the vessel. 

The sample volume is then moved electronically across the vessel in 
discrete steps, with several heart cycles of blood velocity informa- 
tion recorded at each point. The EKG must be recorded simultaneously 
so that velocity measurements can be synchronized in time. Generally 
the quantity of data is such that computer processing of the velocity 
waveforms is required. The requirement that the Doppler angle be 
substantially less than 90® prohibits recording a velocity profile 
across a true diameter; however, for most measurements, the flow 
profiles can be assumed to be constant over the axial length in 
question. If a single gate PUDVM is used, the flow characteristics 
are also assumed to remain constant over the time required. 

In addition to the limitations implied in the above mentioned 
assumptions, use of the diameter profile integration technique is 
subject to other errors. These errors may be divided into two cate- 
gories: 1) experimental errors in applying the technique, and 

2) resolution errors in the measurement of local velocity and vessel 
diameter with the PUDVM. The first category involves mostly experi- 
mental skill and technique while the second category involves inherent 
systematic errors. 
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Certain assumptions are involved in calculating flow from the 
measured profile. First of all, the vessel geometry is assumed to 
by cylindrical over the measurement region. Secondly, the blood 
velocity field should have a degree of radial symmetry about the vessel 
axis. More specifically, the measured velocity profile should be 
valid for rotation about the vessel axis by i 90®. Thirdly, the 
blood velocity vector is assumed to be parallel to the vessel axis 
at all times during the cardiac cycle, maintaining a constant Doppler 
angle to the sound beam. These latter two assumptions imply that 
the profile integration technique is not suitable for highly skewed 
velocity fields. See Daigle (1974) for details of profile integration 
method. 

For transcutaneous measurements, the Doppler angle is usually 
set by two steps. The transducer is first adjusted for a null signal 
which orients the sound beam at 90® to the vessel axis, "he Doppler 
angle is then set by rotating the transducer through a known angle 
in a plane which contains the vessel axis. This two step procedure 
requires a complicated transducer holder for measuring the angles 
involved and is time consuming. Also, refraction effects as the 
sound beam enters the tissue can introduce considerable error with 
this approach. 

3. Profile Truncation and Integration 

We know from theory and experiment that a velocity profile 
recorded with a PUDVM may exhibit tailing at the near and far walls 
of the vessel due to boundary and convolution error. This may cause 
an overestimation of vessel diameter and velocity at the region of 
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the walls. This is the rationale for the max/min slope method and 
reduction of velocity values at the wall described in method (Section 
2). A second flow calculation method is known as profile truncation 
and integration and consists of truncating the measured velocity 
profile to the true vessel diameter (see below). 



profile is integrated as described prev'lously. In these experiments 
the assumption of peak velocity at the centerline was made to permit 
wall location. This method therefore requires an independent, 
accurate measure of vessel diameter. 

4. Diameter Gate Average Velocity Method 

The diameter gate average velocity method utilizes a small 
transducer compared to the size of the vessel to obtain a centerline 
velocity with a narrow gate and a diameter average velocity with a 
wide gate encompassing the full vessel (Daigle, 1974). 
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The mean velocity V Is given by 



where V is the centerline velocity and V is the average velocity 
P ^ 

along a thin sample beam across a diameter. 
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III. VELOCITY SCANS - RESULTS 

1 . Velocity Profile Data 

Nine figures summarize the velocity scan data. Figures 9, 10, 

12, 13, 15, and 16 exhibit scans for various gates for laminar flow 
while Figures 11, 14, and 17 exhibit turbulent profiles. The 
theoretical profile is shown on each figure as a dotted line and the 
measured profiles are indicated with their gate lengths. The trans- 
ducer diameter and Reynold's number are indicated in the upper 
right hand corners of the figures. All data are nondimensional with 
the nondimensional radius on the abcissa and the nondimensional 
velocity on the ordinate. The broadening of the profiles with increasing 
gate length is obvious on all graphs. The boundary error effect 
(Daigle, 1974) is subtle but is more noticeable for gate lengths of 
~ 3R/2 (see Figure 1). These figures should be compared with Figures 
2-8, 2-9, 2-10, and 2-11 (Daigle, 1974). 

2. Velocity Scans - Discussion 

The PUDVM measures centerline velocity accurately for R^i^gnsducer 
< R/2ygggg^ and a gate < R/2ygsger constraint the center- 

line velocity errors are generally less than 435. When the gate is 
increased to R, 3R/2, and 2R the error in centerline velocity is 735, 

1135, and ]B% respectively. Figure 9 shows the velocity profiles 
for transducer R/2 and we observe the profile broadening and Cj^ 
velocity dimunition with increasing gate. The measured profile 
appears somewhat narrower than the actual profile for = R/2. 

For a Reynold's number of 3000 (turbulent), the smaller gates over- 
estimate the velocity by approximately 1 035. The assumed profile 



is quite blunt and decreased bluntness would provide more accurate 
results. For the turbulent case it is difficult to accurately 
choose a profile vis a vis laminar flow (see next section). 

For a transducer diameter R we see a generally decreased accuracy 
in velocity measurements, (Figures 4, 5, and 6) except for the 
turbulent case where the error is - 656. 

For a transducer diameter 2R (full vessel illumination) the 
accuracy of Cj^ velocity nieasurementu drops off still further as 
expected. However, for the turbulent flow case the velocity is 
measured quite accurately for all gates, but again this is for a 
somewhat arbitrary choice of the turbulent velocity profile. 

For a wide gate condition, i.e., Tg ~ 2R the measured center- 
line velocity for different transducer diameters is: 

Velocity 



Dia 

Gate 

Laminar 

Turbulent 

1) 

R/2 

2R 

.82 

1.04 

2) 

R 

2R 

.77 

.97 

3) 

2R 

2R 

.73 

.96 


We are most interested in case 3, the full vessel illumination, wide 
gate condition. Here we find that we overestimate the mean velocity 
by nearly 50%. This result can be attributed to the use of circular 
transducers with overestimation of central core velocities. This 
result is in agreement with all wide gate measurements we have made 
to date using circular transducers. 

Calculation of Turbulent Profile: Reynolds Number = 3000 

From Schlicting (1974 Ed. p. 563) we have the following profile 
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FIGURE 16 
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given for turbulent flow, Re - 3000: 

u(r) = VoO - (1) 

Knowing the flowrate in our experiment (used to determine Re = 3000) 
we can integrate (1) to relate Q, the flow, to the centerline velocity 
for the turbulent profile, 

q = 2uV^ (1 - rdr 

integrating this expression 

Q = ,3956(2TrVQR^) 

therefore 

y = V = ^ — 7- 

0 ,7902ttR^ 

For a flow Q of 18.12 corresponding to a Reynold's number of 3000, 
the centerline velocity is 49.18 cm/sec. However, unlike the laminar 
flow case where we have substantial theoretical and experimental 
background to substantiate the parabolic velocity profile, in tur- 
bulent flow there is more room for error in choosing the profile 
shape as above. Nevertheless, the profile we have chosen should be 
close to the real profile. 
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4. Volume Flow Comparison 

The primary purpose of these experiments was to determine 
the performance accuracy of the PUDVM as a flowmeter. To carry 
out this objective, we varied transducer diameter, gate length, 
and Reynold's number and calculated volume flow by 1) profile 
integration where the wall locations are calculated by the mini/ 
maxi slope method, 2) profile integration truncating the velocity 
at the walls where wall locations are set by knowing the true 
diameter in one case and by measuring the half power diameter 
in the second case, and 3) the small transducer-wide gate method. 

The results of these experiments are summarized in Table I where 
all data has been normalized to actual measured volume flow {by 
rotimeter and volume collection). 

For transducer R/2 increasing the gate from R/2 to 2R produces 
increases in calculated volume flow as expected, but the noteworthy 
result is that the volume flow estimate is most accurate using the 
truncation method regardless of the comparable gate length. This 
result is consistent also for transducers R and 2R. It is also inter- 
esting to note that volume flow can be more accurately measured for 
turbulent flow, where presumably the flattened profile helps improve 
the accuracy by eliminating steep velocity gradients across the 
full vessel lumen. To conclude, tolerable calculations of volume 
flow for all Reynold's numbers were obtained for transducer R/2 
with gates R/2, R; transducer R with gates R/2, R, 3R/2; and trans- 
ducer 2R with gates R, 3R/2, 2R. The R/2 gate for transducer 2R 
could not be obtained due to high initial velocity values. Compare 
these results with data from Table 2-2 (Daigle, 1974), 



VOLUME FLOW COMPARISONS I 






-32- 


5. DIAMETER - Gate Average Velocity 


From Daigle (1974) 

V - 1/2 Vp 

V = 2/3 V 
c ' p 

using data obtained with transducer diameter R/2 
Table I; 

- .82 for gate 2R 
Vp = 1.0 for gate R/2 
therefore calculating V 

■y = V V _ . ...82 _ jQ 

V p C ./u 


V = .70 

Note again this is a 4056 overestimate of V and thus volume flow. Although 
the transducer is slightly larger than R/4, this overestimation is not 
a strong recommendation for the diameter - gate average velocity method. 
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Wide Gate (2R) Centerline Velocities 


Transducer 

RE 

It^M 

^actual 

% overest 

R/2 

1000 

.82 

.50 

44 


1500 

,82 

.50 

44 


3000 

1.04 

.81 

28 

R 

1000 

.75 

.50 

50 


1500 

,78 

.50 

56 


3000 

.97 

.81 

19 

2R 

1000 

.72 

.50 

44 


1500 

.74 

.50 

48 


3000 

.96 

.81 

18 


We can conclude that a rectangular crystal transducer is required for 
accurate wide gate mean velocity measurements in contrast to the circular 
transducer whose results are shown here. 
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IV. CURRENT WORK 

1 . Ultrasonic Transducer Power Emission 

A Cahn Electrobalance model 4600 special has been adapted to 
measure energy emitted by ultrasonic transducers used in this research. 
The purpose of these measurements will permit us to quantitate 
exposure levels in the event the transducers may be applied to human 
patients. 

2. The electronic specifications for the pulsed Doppler electronics 
employing the wide gate method are being set. Fabrication will 
occur following specification of emitted power and wide band signal 
processing. 

3. In vivo flow calibrations are being performed to establish the accuracy 
of the wide gate method for blood flow measurement. 
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